Z00lll ° z ^ 

5 'KM / 


Complex Correlation Kohn-T Method of Calculating Total and Elastic Cross Sections. I: fZ %/Z 

Electron- Hydrogen Elastic Scattering 


A. K. Bhatia and A. Temkin 
Laboratory for Astronomy and Solar Physics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 


Abstract 

We report on the first part of a study of electron- hydrogen scattering, using a method [Bhatia, 
Schneider, and Temkin, Phys. Rev. Lett. 70, 1936 (1993)], which allows for the ab initio calculation 
of total and elastic cross sections at higher energies. In its general form the method uses complex 
“radial” correlation functions, in a (Kohn) T-matrix formalism. The titled method, abbreviated 
CCKT, is reviewed, in the context of electron-hydrogen scattering, including the derivation of the 
equation for the (complex) scattering function, and the extraction of the scattering information 
from the latter. The calculation reported here is restricted to S-waves in the elastic region, where 
the correlation functions can be taken, without loss of generality, to be real. Phase shifts are 
calculated using Hylleraas-type correlation functions with up to 95 terms. Results are rigorous 
lower bounds; they are in general agreement with those of Schwartz [Phys. Rev. 124, 1468 (1961)], 
but they are more accurate and outside his error bounds at a couple of energies. 
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I. INTRODUCTION 

fn a previous paper [l], we developed a method for calculating elastic and total cross sections 
for electron-atom scattering. Called “complex correlation Kohn T’\ and abbreviated (CCKT) the 
method is intended for non-elastic scattering in low and medium energy range, by which we mean the 
energy domain where all the inelastic channels are open plus the continuum of ionization channels 
(up to a total available energy E, which is, say, less than five times the continuum threshold). 
Clearlv no method can specify all such open channels individually, and usual approximations such 
as the Born and eikonal methods [2] are only reliable at high energies. Because our method is 
intended for the low continuum, it is a partial wave method where the partial wave expansion can 
still be expected to converge reasonably effectively. 

The original formulation [l] was applied to an approximation of electron- hydrogen (e-H) scattering 
called the spherically symmetric model [3], The S-wave part of that model reduces to the Temkin- 
Poet model [4,5], and it was only for the latter that the actual calculation was carried out. (The T-P 
model was later generalized in a different and more incisive way in what was called the “generalized 
exchange approximation” [6] ). ° 

In this paper we deal with the real (i.e. physical ) e-H problem. The method is reviewed in 
Sec II. Specifically a more detailed derivation of the (complex) optical potential and scattering 
equations is given. As opposed to Ref. [l], the emphasis here is on the numerical solution of 
the scattering equation by decomposing it into its coupled real and imaginary parts and from the 
solution extracting the (non-unitary) T-matrix. 

An actual calculation in its full generality, however, is still a very extensive undertaking. We shall 
therefore confine ourselves in this initial study, to elastic S-wave scattering. Here our aim is to 
check and increase the accuracy of the classic calculation of Schwartz [7]. Results are presented 
in Sec. III. Suffice it here to say that our results are rigorous lower bounds on the phase shifts. 
Thus where they exceed Schwartz’ results, which they do at two energies, they are necessarily more 
accurate. In general our results are within the stated uncertainty of his remarkable calculation 
(considering when it was done). Here, given the convergence of our results, we would only claim 
greater accuracy. 

II. THE COMPLEX CORRELATION KOHN T (CCKT) METHOD 
A. General 

Confined to the e-H partial wave (denoted by L) problem, the total spatial function is written: 

^^(rirTz) = ^^Y LQ {r\)4> lo {r 2 )Yoo{r\) ± (1 2) + ^(n, r 3 ) (I) 

The superscript (±) above refers to singlet (upper sign) or triplet (lower sign) scattering respec- 
tively. [That superscript will be dropped hereinafter, except where it is necessary.] Beyond the 
terms containing ui explicitly (those are the terms giving rise to the exchange approximation), the 
function is the correlation function, which is our concern here. For arbitrary L this function is 
most efficiently written in terms of symmetric Euler angles [8] : 

H 2T [/r ,£ (n, r 2 , r l2 )Vl' c ($, 4> , 0) + fZ’~'(r 2 . r,. r, 2 )Z>£-' (0 , 0, 0)] (2) 

The c functions (c=+l,-l) are exchange adapted combinations of Wigner functions. (The V 
functions, in general have been called rotational harmonics [9].) The / ’s above are generalized 
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“radial” functions, which depend on the 3 residual coordinates that, are required (beyond the Euler 
angles) to define the two vectors iq and rj (Further technical details of the symmetry and other 
properties of the symmetric Euler angle description are given in [$]). What distinguishes this CCKT 
approach is the fact that the “radial” functions are taken to be complex. That is different from 
ordinary variational methods, wherein the radial functions are taken to be real. In the latter case, 
when one forms the variational functional to be varied [cf. Eq.(5j below], one uses the function 
) on the left. That is so because the / functions would be real, so that only the V functions in 
(2) would be complex conjugated. When the radial functions are complex, however, then - in order 
to have a variational principle for the T matrix - only the P but not the radial functions are to be 
complex conjugated in the left-hand function. That is indicated by using the tilde rather than star 
for that function, which is used in the functional I L of the variational principle (below). Explicitly 


q> L = q> L - 


(3) 


where 

$L = y] yi [/£ ,e ( r l. r 2 » r ! 2 ) 4>, $))* + /£’~ e ( r 2 . ri, r 12 ) (x>£ ‘(0,0,^)) ] (4) 

£ K 

. The Kohn variational principle [10] for Ti reads : 


s[i l -(-1) l T l ]=0 (5) 

where II is the (Kohn) functional 

I L =[ [*l{H - E)9 l ] d ri dr 2 (6) 

Carrying out the variation leads to scattering equation for u^^r) (Rydberg units throughout): 


-dl _ + Vj ± v„ + vw - 

V d and V ex are the well known direct and exchange potentials of the “exchange approximation” [11] 
Those potentials (the latter being non-local) are real, however the (remaining) optical potential 
acting on U£,(r) is complex : 


,(±) - 


= 0 


(7) 


= r , (8) 

Thus, ui(r) is nontrivially complex, which leads to a nonunitary Tl matrix, which is derivable 
from the asymptotic form of ui 


lim 

r— ►oo 


UL{r) 


i L $in{kr - ttZ,/2) + j^ikr 


(9) 


In defining the optical potential we have used the Feshbach approach [12], involving projection 
operators P and Q, which for the hydrogenic (i.e. one-electron) target can be written expicitly [13] 
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P = Pi + P 2 - PlP'l 


( 10 ) 


where the spatial projectors are : 


Q =1- P 


(U) 


P, = e r 'Y 00 (ri) >< e r 'Vbo(f.) (12) 

Note, Pi and P 2 commute and are each idempotent, hence the complete P and Q operators are 
idempotent (P 2 =P; Q 2 —Q) and orthogonal (PQ= 0). 

To complete this review of the CCKT method, the optical potental is expanded in terms of the 
eigenspectrum of the QHQ problem: 
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< > 

< $lQ$l > 


= 0 


(13) 


This leads to complex radial eigenfunctions and complex eigenvalues £ s . Insertng a complete 
set of these functions (understood to be approximated by a discrete set arising from (13) using an 
ansatz with a finite number of variational parameters) allows the optical potential to be expanded 


5 tr — C S 

To repeat the main point of the CCKT method: because the radial functions /£’ e (ri, r 2 , r i2 ) in (2) 
are complex, the resulting Ti matrix (which in this method is actually a number) will be complex 
and nonunitary. This means that the elastic cross section, calculated from 


a L = 4w(2L+l)\T L \ 2 (15) 

will be different from the total cross section 

47 r 

a L = —(2 L + 1 )Im{T L ) (16) 

The effect of a nonunitary Tl, is of course well known from the literature (cf., for example Ref.[14]). 
What is new about the CCKT is that it provides an ab initio method for calculating X^,, rather than 
parametrizing it, as was done for example in the ’’clouded crystal ball” model of neutron-nucleus 
scattering and reactions [15]. 

B. Some Details of the CCKT Method 

The Kohn variational principle is generally applied in such a way, as in Ref.[l], that the function 
Ui.fr) of Eq.(7) is parametrized by an expansion in terms of known functions, and the functional 
in Eq.(fi) is analytically evaluated; it is a bilinear expression in terms of the (linear) expansion 
coefficients. The (Kohn) principle then reduces to differentiation with respect to these parameters, 
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which in turn leads to a set of linear equations, the last being an inhomogeneous equation from the 
differentiation with respect to Tl- The solution of the set of equations then provides values for all 
the (linear) constants plus the value of Tl- The latter can be improved by subtracting the value 
of the functional itself, which in general will not be zero when the calculation is performed in this 
way. In that case an improved Tl can be obtained by subtraction: 

^improved -t l - r L ( 17 ) 

When, however, the calculation is carried out (as it is here) by solving directly for u L (r), then this 
is equivalent to making Il= 0. Thus the Tl that emerges is automatically the improved (i.e. the 
second order) result. What we shall now do is to outline briefly how the calculation is performed 
when the potentials and hence the solutions are complex. Specifically, the radial equation can be 
written: 



(18) 


where 

4« = _^ + Mi+ii + ^±v; I -t 2 (W) 

The V d is the (local) Hartree potential [ V d {v) = -2 exp(-2r) (1 + 1/r)] and V ex is the well-known 
(non-local) static exchange potential [11]. The remaining (non-local) potential is what is gererally 
called the optical potential, Eq.(5). 

We now give a few details of how one solves the radial equation (18) in the general case. Dropping 
sub (super) scripts, we write the solution of (4) 


u(r) = uji(r) + iui(r) (20) 

where u^(r) and u/(r) are real functions which, substituted in (18), lead to a coupled set of (real) 
differential equations. 

(£ + vW)u R = (21) 

(£ + V^)ui = -V^ur (22) 

Vi? and V^Jp above are the real and imaginary parts of the optical potential, respectively. (Both 
are real, non-local potentials). The coupled Eqs.(2l, 22) have two sets of solutions (i.e. both regular 
at the origin), labeled u^(r), uj°( r) [i = 1,2]. From these two solutions one determines two real 
constants C\ and C 7 such that the real and imaginary parts of the asymptotic form (9) is satisfied, 
that the equation for L = 0 reduces to 


lim 

r-*oc 


C' n ^ 

CiUfl 


+ C 2^ 


( 2 ) 

n 


sinkr 

k 


+ TReoskr - Tjsinkr 


lim + C -2 u /* ^ = T jC.oskr + TRsinkr 

r-*<Xi l J 


(23) 


(24) 



The two linearly independent solutions will necessarily have the asymptotic form 


, l l^ u R = A R*‘n(kr+T}%) 

(25) 

r !i& “/ 1 = A l i ] sin(kr + r?} 0 ) 

(26) 

in which one of the normalization constants, say ,4^\ is arbitrary. 
,4 (l )=A^/V-4^\ which allows (23, 24) to be rewritten 

This is tantamount to letting 

CiCOST]^ + C2COST}^ = i - Tj 

K 

(27) 

Cisinrj^ C 2 sinr$ = Tr 

(28) 

CiA^cosrfp + C 2 A^ cosrf-p = Tr 

(29) 


CiA^sin^ + C 2 A^ sinr)[ 2 ^ = 7/ (30) 

The T]'s and A’s are extracted from the solutions of (21, 22); thus (27-30) is a set four equations 
for four real constants : C\ , C 2 and in particular the real and imaginary parts of T: Tr and Tj 
From the latter one calculates the elastic and total cross sections, Eqs. (15), (16). 

III. CALCULATION 

As stated in the INTRODUCTION we calculate here only S- wave (i.e. L=0) elastic scattering, 
specifically phase shifts. This means, since V L = constant, that the correlation function $£, is only 
a function of the “radial” coordinates. Here $£ is taken of Hylleraas form 

;v„ 

*So = e- 7r ‘- 5r2 £ Ci mn r[r?r n l2 ±(1«2) (31) 

Iran 

where the sum includes all triples such that l ttx- j- ti—& and 0, 1, 2, ....... 8 . The total number 

of terms N„ depends on spin and whether 7 =5 or not. The values of N„ as a function of /, m, n is 
given in the following simple array: 

w-* 01234 5 6 7 8 

'S: N w (-r~S) 1 3 7 13 22 34 50 70 95 

3 S: N w ( 7 # 5) 1 4 10 20 35 56 84 120 165 

To summarize the calculation: the QHQ problem is solved (for a given 7 and S and N w ). The 

result is a set of eigenvalues £ s (s= 1,2, .N w ) and associated eigenfunctions From them 

(the terms in) the optical potential, Eq.(14) is constructed, and .the integro-diflerential Eq.(7) is 
solved noniteratively. Because all quantities are real, the solution is unique (up to an arbitrary 
normalization) with asymptotic form 

G 



(32) 


lim a sin(kr + q^ ±] ) 

From u (±) (r) and as derivative the phase shift q is readily extracted. 

Examples demonstrating the convergence of q for A:=0.S as a function of iV w are given in Table I. By 
virtue of the fact that q 1 s are rigorous lower bounds on the phase shift [13], the convergence then 
becomes a powerful indication of the accuracy of the result. Phase shifts as a function of lc are given 
in Tables II and III. To repeat: they are rigorous lower bounds in all cases. They are compared 
to the results of Schwartz [7]. Our l S phase shifts are seen to exceed Schwartz' (including his 
estimated error) at Ar=0.3. 0.4. Otherwise our results for both l S and 3 5 are within his estimated 
uncertainty, which we find to be quite impressive considering that his calculation was carried out 
over 40 years ago! Nevertheless we believe that our present calculation is more accurate. The 
difficulty in making a secure extrapolation is due to the difficulty of estimating quantitatively the 
effects of polarization and other long-range potentials at non-zero k. (In principle, one can optimally, 
rigorously include polarization by use of a polarization pseudostate [16]; however in practice those 
corrections are small, yet require a much more extensive calculation than is worthwhile here for 
our present purpose.) The convergence of the results suggests that they are accurate to one or 
two units in the fourth place after the decimal point. (We have given five digits because to that 
accuracy they are rigorous lower bounds.) 

Scattering lengths are not included in our calculation, because at fc=0 the Kohn variational prin- 
ciple, which applies to Schwartz’ calculation, is well known to provide rigorous (in this case upper) 
bounds [17]. It is worth recalling that one of us [18] had earlier deduced that long range polarization 
would have decisive quantitative effect on the scattering length and had derived a formula 

a-.(«)-«(i--^ + 0(i)) (33) 

in which a(R) represents the value of scattering length including only the region of configura- 
tion space r < R, and the remaining terms in (33) give the effect of the long range polarization 
(a=polarizabiIity of the hydrogen atom= 4.5). Using that formula led to the first correct estimate 
of the true scattering length [18]. Thus when Schwartz applied his calculation at k= 0, he found it 
necessary to include the equivalent of long range polarization terms, as predicted [18], in addition 
to the Hylleraas terms in order to get adequate convergence. His results are: a( 1 5)= 5.965 ± 0.003, 
a( 3 S)= 1.7686 ± 0.0002. For these reasons we believe his scattering lengths are sufficiently secure 
and accurate as to not require further calculation. 

Numerical results were obtained with Cray Y-MP computer of the NASA Center for Computation 
Science. 
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TABLE I. Convergence’ of S-wave phase shift as a function of N(~) for A: = 0 . 8 



.... • f 

5 


3 5 


N(~0 

7 = & 

ijl+J 

N(w) 

7 

6 

7 ( ] 

E.A* 


0.65127 

EA h ~ 



1.61729 

3(1) 

1.23 

0.79028 





7(2) 

0.68 

0.87536 

10 ( 2 ) 

0.65 

0.55 

1.63S45 

13(3) 

0.68 

0.87684 

20(3) 

0.84 

0.70 

1.64205 

22(4) 

0.69 

0.88362 

35(4) 

0.85 

0.50 

1.64294 

34(5) 

0.61 

0.88584 

56(5) 

0.85 

0.50 

1.64344 

50(6) 

0.54 

0.88616 

84(6) 

0.85 

0.50 

1.64379 

70(7) 

0.66 

0.88687 





95(8) 

0.72 

0.88718 






° The phase shifts (in radians) have been optimized with respect to 7 and 5 for each N(u>). 
b E. A. are the well-known exchange approximation phase shifts [ 11 ]; it corresponds to no 
correlation terms [N(u>)=0 — > V op =0 in Eq.(7)]. 
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TABLE II. Phase shifts of l S for various k for N=95 


k 

7 = 6 

V 

VSchu/artz 

0.1 

0.53 

2.55358 

2.553(1) 

0.2 

0.60 

2.06678 

2.0673(9) 

0.3 

0.72 

1.69816 

1.6964(5) 

0.4 

0.70 

1.41540 

1.4146(4) 

0.5 

0.70 

1.20094 

1.202(1) 

0.6 

0.70 

1.04083 

1.041(1) 

0.7 

0.71 

0.93111 

0.930(1) 

0.8 

0.72 

0.88718 

0.886(1) 
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TABLE III. Phase shifts of 3 S for various k for 


k 7 S t] 


'HSchwartz 


0.1 

0.85 

0.12 

0.2 

0.84 

0.39 

0.3 

0.90 

0.37 

0.4 

0.91 

0.38 

0.5 

0 . 8 S 

0.44 

0.6 

0.77 

0.52 

0.7 

0.70 

0.55 

0.8 

0 .S 6 

0.50 


2.93853 2 . 9388 ( 4 ) 

2.71741 2 . 7171 ( 5 ) 

2.49975 2 . 4996 ( 8 ) 

2.29408 2 . 2938 ( 4 ) 

2.10454 2 . 1046 ( 4 ) 

1.93272 1 . 9329 ( 8 ) 

1.77950 1 . 7797 ( 6 ) 

1.64379 1 . 643 ( 3 ) 







